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ABSTRACT 

PSR J1846—0258 is an object which straddles the boundary between magnetars and rotation pow¬ 
ered pulsars. Though behaving for many years as a rotation-powered pulsar, in 2006, it exhibited 
distinctly magnetar-like behavior - emitting several short hard X-ray bursts, and a flux increase. 

Here we report on 7 years of post-outburst timing observations of PSR J1846—0258 using the Rossi 
X-ray Timing Explorer and the Swift X-ray Telescope. We measure the braking index over the post- 
magnetar outburst period to be n = 2.19 ± 0.03. This represents a change of An = —0.46 ± 0.03 or 
a 14.5 cj difference from the pre-outburst braking index of n = 2.65 ± 0.01, which itself was measured 
over a span of 6.5 yr. So large and long-lived a change to a pulsar braking index is unprecedented and 
poses a significant challenge to models of pulsar spin-down. 


1. INTRODUCTION 

Many of the quoted properties of pulsars, such as the 
surface magnetic field, the characteristic age, and the 
spin-down luminosity, are based on the assumption that 
pulsars are well modeled as a magnetic dipole in a vac¬ 
uum. One of the ways we have to test the validity of this 
assumption, and by doing so probe the emission mecha¬ 
nisms of pulsars, is by measuring the change of a pulsar’s 
spin-down rate over time. This is expected to behave fol¬ 
lowing a power law, 

v = —Ku n , (1) 

where v is the spin frequency of the pulsar and n is re¬ 
ferred to as the ‘braking index’. In the canonical case of a 
rotating magnetic dipole in a vacuum, the braking index 
is expected to be 3 (e.g. Manchester & Taylor 1977). 

Observationally, the braking index is measured by 
means of observing a gradual change in the frequency 
derivative, and expressing 


Measuring braking indices for pulsars has proven to be 
difficult, as young pulsars typically exhibit large amounts 
of timing noise which can contaminate measurements of 
i) (e.g. Hobbs et al. 2010). As yet, only eight pulsars have 
measured braking indices (see Roy et al. 2012; Espinoza 
et al. 2011, and references therein) which range from 0.9 
(Espinoza et al. 2011) to 2.91 (Weltevrede et al. 2011). A 
braking index of —1.5 was reported for PSR J0537—6910 
but the timing behavior is complex, and dominated by 
glitches as discussed in Middleditch et al. (2006). 

PSR J1846—0258 is a ~800 year old pulsar located 
in the Kesteven 75 supernova remnant (Gotthelf et al. 
2000). It has a rotation period of ~327 ms and is one 
of the youngest known pulsars. For the majority of 
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its observed lifetime, PSR J1846—0258 behaved as if 
it were a typical rotation-powered pulsar, with its X- 
ray emission being much less than the luminosity ex¬ 
plainable by its spin-down power. Curiously, however it 
has no detectable radio emission (Archibald et al. 2008). 
PSR J1846—0258 is also one of the eight pulsars with a 
measured braking index, observed to be 2.65±0.01 from 
2000 to 2006 (Livingstone et al. 2006). 

In 2006, PSR J1846—0258 underwent a rare event - 
its pulsed X-ray flux increased dramatically, it had a 
large glitch, and emitted several magnetar-like bursts 
(Gavriil et al. 2008; Kumar & Safi-Harb 2008; Kuiper 
& Hermsen 2009). PSR J1846—0258 remains the only 
seemingly rotation-powered pulsar to display such dis¬ 
tinctly magnetar-like behavior, making it an interesting 
transition object between the two classes. 

After this magnetar-like outburst, PSR J1846—0258 
went back to manifesting itself as a rotation-powered 
pulsar (Livingstone et al. 2011). However, after tim¬ 
ing the source for more than two years post-outburst, 
Livingstone et al. (2011) measured a braking index of 
n = 2.16 ± 0.13 during this period, a value inconsis¬ 
tent with the braking index measured prior to the out¬ 
burst. While the braking index is expected to change on 
a timescale of thousands of years (see e.g. Contopoulos & 
Spitkovsky 2006; Gourgouliatos & Cumming 2015), such 
a sudden change is unexpected in the standard models. 

Here we report a further five years of X-ray timing 
observations of PSR J1846—0258, for a total of seven 
years after the magnetar-like outburst. We show that 
the braking index is consistent with the post-outburst 
measurement of Livingstone et al. (2011), and inconsis¬ 
tent with that prior to the outburst. This indicates that 
the 2006 magnetar-like outburst resulted in a persistent 
change in the braking index in the source. 

2. OBSERVATIONS AND ANALYSIS 
2.1. RXTE 

In this work, we analyze observations of 
PSR J1846—0258 from the Proportional Counting 
Array (PCA) aboard the Rossi X-ray Timing Explorer 
(RXTE) from January 2008 until the decommissioning 
of RXTE in 2011 December. The PCA consists of five 
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collimated xenon/methane multianode proportional 
Counter units (PCUs) which are sensitive to photons 
in the 2-60 keV range (Jahoda et al. 1996, 2006). The 
PCA was operated in “Good Xenon” mode, which 
provides l-/is resolution for photon arrival times. 

Observations were obtained from the HEASARC 
archive and barycentered to the location 
of PSR J1846—0258, RA = 18 /l 46 m 24.94 s , 

DEC = —02°58'30.1" (Helfand et al. 2003) using 
the barycorr tool in HEASOFT u 6.16. Observations 
were filtered to remove non-astrophysical events using 
xtefilt. In order to maximize the signal-to-noise ratios 
of pulse profiles so as to minimize uncertainties on 
resulting pulse times-of-arrival (TOAs) (see §3), we used 
events from all layers of the then-operational PCUs. 

In total 363 RXTE observations providing ~ 0.9 Ms of 
exposure time were analysed in this work spanning Jan¬ 
uary 2008 to December 2011. Observations taken within 
2 days of each other were merged, resulting in 177 TOAs 
for a typical exposure time of 5-ks per TOA. 

2.2. Swift XRT 

We began observing PSR J1846—0258 with the Swift 
X-ray Telescope (XRT) on 2011 July 25 as part of a cam¬ 
paign to monitor several magnetars (see e.g. Archibald 
et al. 2013; Scholz et al. 2014; Archibald et al. 2015b). 
The Swift XRT is a Wolter-I telescope with a e2v CCD22 
detector, sensitive in the 0.3 — 10 keV range. The XRT 
was operated in Windowed-Timing (WT) mode for all 
observations. This gave a time resolution of 1.76 ms. 

Level 1 data products were obtained from the 
HEASARC Swift archive, reduced using the xrtpipeline 
standard reduction script, and barycentered to the loca¬ 
tion of PSR J1846—0258, using HEASOFT u6.16. In¬ 
dividual exposure maps, spectra, and ancillary response 
files were created for each orbit and then summed. We 
selected only Grade 0 events for spectral fitting as higher 
Grade events are more likely to be caused by background 
events (Burrows et al. 2005). To maximize the signal-to- 
noise ratios of pulse profiles so as to minimize uncer¬ 
tainties on resulting pulse times-of-arrival (see §3), only 
photons from 2.7-10 keV were used. 

To investigate the flux and spectral evolution of 
PSR J1846—0258, a circular region having a 10-pixel ra¬ 
dius centered on the source was extracted. As well, an 
annulus of inner radius 75 pixels and outer radius 125 
pixels centered on the source was used to extract back¬ 
ground events. 

In total 66 XRT observations totaling 541 ks of expo¬ 
sure time were analyzed in this work. Observations taken 
less than 5 days apart were grouped to extract a single 
TOA yielding 47 TOAs, with a typical exposure time of 
10 ks per TOA. 


3. TIMING ANALYSIS 
3.1. Phase-Coherent Timing Analysis 

TOAs for all RXTE and Swift observations were ex¬ 
tracted using a Maximum Likelihood (ML) method as 
described in Livingstone et al. (2009) and Scholz et al. 
(2012). The ML method compares a continuous model 
of the pulse profile to the photon arrival times obtained 
from a single observation. In order to create the con¬ 
tinuous model of the pulse profile, first we create a high 
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Fig. 1. - Timing residuals of PSR J1846—0258 from MJD 54492- 
56880 (post-outburst) for the solutions presented in Table 1. The 
top panel shows the residuals of Solution 1. The bottom panel 
shows the residuals of Solution 2. The vertical dashed line indicates 
where there is a phase ambiguity; see §3 for details. Note that 
Solution 1 is fitted only to data before the phase ambiguity, and 
Solution 2 only to those after. 

signal-to-noise template profile by folding many obser¬ 
vations together using a whitened timing solution. For 
RXTE, the template was derived from folding all pre¬ 
outburst observations, and for Swift, using all the ob¬ 
servations. Separate templates were used for the RXTE 
and Swift observations to account for differences in the 
responses of the telescopes. In both cases, a continu¬ 
ous model of the profile was created by fitting the high 
signal-to-noise template with a Fourier model using the 
first two harmonics. Two harmonics were chosen to op¬ 
timally describe the pulse shape, as determined by the 
h-test (de Jager et al. 1989; Archibald et al. 2015a). 

These TOAs were fitted to a timing model in which 
the phase as a function of time t can be described by a 
Taylor expansion: 

= (fo+Vo{t-to) + 7}Vo{t-t 0 ) 2 +-P 0 (t-t 0 ) 3 +- ■ ■ (3) 

where vq is the rotational frequency of the pulsar at time 
to- This was done using the TEMP02 (Hobbs et al. 2006) 
pulsar timing software package. 

In Figure 1 we show the timing residuals in the range 
MJD 54492 to 56880, the period after the magnetar-likc 
outburst and glitch recovery have relaxed. For details 
about the glitch and the glitch recovery, see Kuiper & 
Hermsen (2009) and Livingstone et al. (2010). 

Finding a single phase-coherent solution over the entire 
seven-year post-outburst data set is not possible due to a 
phase ambiguity during the Sun constraint period from 
MJD 56246 to 56338. This is indicated in Figure 1 by 
a dashed vertical line. We were able to find two phase- 
coherent solutions, one before this Sun constraint and 
one after. The two timing solutions are presented in 
Table 1. 
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TABLE 1 

Phase-Coherent Timing parameters for PSR J1846—0258. 


First Phase-coherent Solution 

Dates (MJD) 

54492.0-56246.7 

Dates 

2008 Jan 27 - 2012 Nov 15 

Epoch (MJD) 

55369.00000 

* W 1 ) 

3.059 040 903(4) 

* W 2 ) 

-6.65131(1) x lO” 11 

V (s -3 ) 

2.937(8) X 10“ 21 

rms residual (ms) 

304.4 

rms residual (phase) 

0.931 

Braking index, n 

2.031(6) 

Second Phase-coherent Solution 

Dates (MJD) 

56338.7-56964.20 

Dates 

2013 Feb 15 - 2014 Nov 03 

Epoch (MJD) 

56651.00000 

* W 1 ) 

3.051 693 972(3) 

V (s~ 2 ) 

-6.61349(2) x 10“ u 

V (s- 3 ) 

3.30(4) x 10~ 21 

rms residual (ms) 

32.4 

rms residual (phase) 

0.099 

Braking index, n 

2.30(3) 


Note: Figures in parentheses are the nominal 1<j tempo2 
uncertainties in the least-significant digits quoted. 


This loss of phase coherence could be due either to 
a glitch, or to timing noise. Fitting for a glitch during 
the Sun-constraint period using a timing solution up to 
i'0 yields Av/v = 5.7 ± 0.5 x 10 -8 and A!>/;> = —2.5 ± 
0.4 x ICE 4 over the Swift campaign. Fitting using both 
the RXTE and Swift data sets gives glitch parameters 
ranging from Av/v of —9 x 10“ 8 to 1.7 x 10~'. We note 
that these values vary based on the time-span fit and the 
number of frequency derivatives used in the fit. Finally, 
we note that fitting a continuous solution over the gap 
yields comparable residuals to the glitch fits. Thus we do 
not need to invoke a sudden glitch to explain the timing 
behavior of PSR J1846—0258 at this epoch. 

3.2. Partial Phase-Coherent Timing Analysis 

Measurements of v can be susceptible to contamination 
from timing noise (e.g. Hobbs et al. 2010). To mitigate 
this effect, we fit small segments of data to make local 
measurements. For all methods presented below, rela¬ 
tive pulse numbers were fixed to those given by the fully 
phase-coherent timing solution. No timing solution was 
fit overlapping the phase ambiguity. 

For each small segment of data, using the established 
pulse numbers, TOAs were fit to a timing solution con¬ 
sisting of only v and v. The time spans were determined 
by allowing a maximum xt of ~1 and the condition 
that there was no apparent-by-eye red-noise signal in the 
residuals. When this condition was met, we moved over 
by half the number of TOAs in that solution, and fit 
again until the criteria were met. We did not allow a 
solution to span over a Sun-constraint period. 

In Figure 2, we show these measurements of 0 over the 
data set. The top panel shows z> over time. The mid¬ 
dle panel shows v over time subtracting a constant slope 
consisting of the pre-outburst braking index, 2.65±0.01 
(Livingstone et al. 2011). Note the clear linear trend in 
the middle panel indicating that the pre-outburst brak¬ 
ing index does not describe the data well. 

We fit a slope to the post-outburst v in order to obtain 
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Fig. 2.— v measurements for PSR .11846—0258 from MJD 
54492-56880. The top panel shows the measured is. The solid 
black line shows the pre-outburst v of 3.88 X 10 -21 s '’’.The mid¬ 
dle panel shows the same data subtracting the pre-outburst v. 
The black line in this panel shows the difference between the 
pre-outburst measurement, and the best-fit post-outburst v of 
3.17 ± 0.05 X 10 -21 s -3 . The gray shaded region shows the 1- 
a bounds on this determined from a bootstrap analysis to the full 
data set as described in the text. The bottom panel shows the ;> 
residuals after subtracting the best-fit slope from above. The ver¬ 
tical dashed line indicates where there is a phase ambiguity; see §3 
for details. 

a measurement of v , and thus a braking index. The 
timing measurements of PSR J1846—0258 have a scatter 
larger than would be suggested by their formal errors, 
therefore we use a bootstrap method. The bootstrap 
method is robust for error estimation when only a small 
number of measurements are available (Efron 1979) and 
the formal uncertainties are thought to not fully describe 
the data. 

For the full post-outburst data set, this yielded a mea¬ 
surement of v — 3.17 ± 0.05 x ICE 21 s -3 corresponding 
to a braking index of n = 2.19 ± 0.03 for the bootstrap 
method. The residuals of this fit can be seen in the bot¬ 
tom panel of Figure 2. 

In order to verify that the phase ambiguity between 
the two timing solutions presented in Table 1 does not 
affect our result, we split the data into the corresponding 
two segments. Fitting from MJD 54492.0-56246.7 gives 
a braking index of n = 2.09±0.05. Fitting MJD 56338.7- 
56880.5 gives n = 2.23 ± 0.07. These two segments gave 
consistent slopes at the 1.4-cr level, and are both incon¬ 
sistent with the pre-outburst braking index. This gives 
us confidence that the measured post-outburst braking 
index ofn = 2.19±0.03 represents a long-lived change 
in the braking index of An = —0.46 ± 0.03, a 14.5 a 
difference. 


Year 

2008 2009 2010 2011 2012 2013 2014 



3.3. Timing Noise 

In order to quantify the effect of timing noise that could 
be contaminating the measurement of the braking index, 
we fit a timing solution consisting of a frequency and 
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Fig. 3.— Timing noise in PSR J1846—0258 over 15 years of X-ray 
timing as described by the A - parameter, see equation 4 . The 
vertical dashed line indicates the epoch of the outburst. 

three frequency derivatives for each year, ending a solu¬ 
tion at times of a glitch, or the start of Sun-constraint. 

Following the method of Livingstone et al. (2011), we 
measure the quantity 


Ay = lOg 


V24 * ) 


( 4 ) 


where t is the length of time over which the solution was 
fit, ~ 2.5 x 10 7 s. This is analogous to the A g parameter 
of Arzoumanian et al. (1994) where As is used as an 
estimation of the contributions of u to the accumulated 
phase deviation of the pulsar. As i) is physically relevant 
in timing measurements of PSR J1846—0258, we use Ay 
as an estimate of the phase contamination from z>' and 
higher order effects. 

In Figure 3, we show Ay over the 15 years of timing 
of this source. While the scatter is high, Az>' shows a 
possible increase for the period following the magnetar- 
like outburst in 2006. 

Before the outburst, eg. from 2000 to 2006, the 
weighted mean was Ay = 0.1 ± 0.2. For the first pe¬ 
riod after the outburst, 2007, Ay = 1.16 ± 0.03, sub¬ 
stantially higher than at any other time. After this, the 
timing noise decreased to a level that is marginally higher 
than the pre-outburst noise, with the weighed mean of 
A y = 0.6 ±0.2 from 2008-2014. Thus, the level of timing 
noise clearly increased following the magnetar-like out¬ 
burst but appears to be relaxing back to the pre-outburst 
level on a time scale of several years. 


4. RADIATIVE PROPERTIES 
4.1. Spectral Analysis 

Swift XRT spectra were extracted from the selected re¬ 
gions using extractor, and fit using XSPEC package ver¬ 
sion 12.8.2 5 . Spectral channels were grouped to 1 count 
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per bin, and fitted using cstat minimization. The spec¬ 
trum was fit with a photoelectrically absorbed power law. 
Photoelectric absorption was modeled using XSPEC tbabs 
with abundances from Wilms et al. (2000), and photo¬ 
electric cross-sections from Verner et al. (1996). Due to 
both the nature of the windowed timing read-out mode 
of the XRT, and the fact that the XRT point spread 
function is comparable to the size of the bright, central 
region of the nebula, we are unable to separate the flux 
coming from the pulsar itself from the bulk of the pulsar 
wind nebula which surrounds it. 

As all of the Swift observations had consistent flux 
and spectral parameters, we co-fit all observations si¬ 
multaneously. This yielded a best-fit model with Nh = 
(4.43 ± 0.05) x 10 22 cm" 2 and P = 1.80 ± 0.02. We note 
that the best-fit power-law index is consistent with that 
of the pulsar wind nebula reported by Kumar & Safi- 
Harb (2008) and Ng et al. (2008), as well as that reported 
in the 20-300 keV range using INTEGRAL (Kuiper & 
Hermsen 2009). 

The absorbed 0.5-10 keV X-ray flux measured over 
the Swift campaign of the combined pulsar and pulsar- 
wind-nebula was (2.04 ± 0.02) x 10” 11 erg cm -2 s _1 . 
In Chandra observations taken 2000 (Ng et al. 2008) 
and 2009 (Livingstone et al. 2011), the absorbed 0.5- 
10 keV flux from the combined pulsar and pulsar- 
wind-nebula were (1.81±0.03)x 10~ n erg cm~ 2 s _1 and 
(1.73±0.07) x 10 -11 erg cm -2 s _1 respectively. While for¬ 
mally, our measured Swift flux and the archival Chandra 
fluxes are inconsistent, the cross-calibration between X- 
ray instruments is only accurate to the ~10% level; see 
Tsujimoto et al. (201l). Therefore we find no evidence of 
a changing flux for the system to the level of instrumental 
uncertainties. 

4.2. Pulse Profile Analysis 

To look for changes in the pulse profile, we folded each 
observation using 16 phase bins. Each profile was com¬ 
pared to the high signal-to-noise-ratio pulse template de¬ 
scribed in §3 by subtracting a fitted DC offset, and using 
a multiplicative scaling factor to minimize the difference 
between the template and scaled profile as determined by 
a x 2 minimization. For both RXTE and Swift , all profiles 
are consistent with the respective telescope’s standard 
template. 

It has been shown previously that PSR J1846—0258 
exhibited no significant change in its X-ray pulse pro¬ 
file during the magnetar event (Kuiper & Hermsen 2009; 
Livingstone et al. 2010). For RXTE we now have compa¬ 
rable data from both before and after the magnetar-like 
outburst. This allowed us to search for long-term lower- 
level changes in the pulse profile. To do so, we combined 
all observations for which we had a valid timing solu¬ 
tion into a two high signal-to-noise-ratio profiles with 64 
phase bins representing the pulse profile before and after 
the magnetar-like outburst. To do this, we fitted each 
year of TOAs to a timing solution, using as many fre¬ 
quency derivatives as necessary to whiten the residuals. 
Each year of data was then folded into a 64-bin pro¬ 
file, and aligned with other years’ profiles using cross¬ 
correlation. This resulted in two high signal-to-noise- 
ratio profiles: the first using 918 ks of exposure time 
from January 2000 to March 2006, and the second using 
871 ks from January 2008 to December 2011. These two 
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Fig. 4.— Normalized RXTE pulse profiles of PSR J1846—0258. 
The solid black profile shows the profile from January 2000 to 
March 2006, just before the 2006 outburst. The red dotted profile 
shows the profile from 2008 January to 2011 December. The bot¬ 
tom panel shows the residual difference between the two profiles. 
The residuals have x^/(d°f) = 0.988/(62) indicating the profiles 
are statistically identical. 

normalized, DC-subtracted, high signal-to-noise profiles 
are shown in Figure 4, as well as the difference between 
them. The residuals have xl/(dof) = 0.988/(62) indi¬ 
cating the profiles are statistically identical. This is con¬ 
sistent with the lack of profile change reported by Kuiper 
& Hermsen (2009) and Livingstone et al. (2011). 

4.3. Burst Search 

All Swift observations were searched for magnetar-likc 
bursts by binning the source region light curves into 0.01- 
s, 0.1-s, and 1.0-s bins. The counts in each bin were 
compared to the mean count rate of its Good Timing 
Interval (GTI), assuming Poisson statistics, similar to 
the methods described by Scholz & Kaspi (2011). We 
found no significant bursts in the Swift observations. 

For the RXTE PCA, due to the background being 
highly variable, each 60-second interval was treated simi¬ 
larly to a Swift GTI. An additional constraint was placed 
on the PCA data that a putative burst must be detected 
in all operational PCUs to be considered real. We find a 
previously unreported burst on MJD 55070, 27 August, 
2009. This burst has a Tgo, the time duration in which 
90% of a burst’s fluence is collected, of 7 ± 1 ms and a 
fluence of 12 ± 3 counts per PCU (24 ± 5 total counts). 
This corresponds to a false alarm probability of ~ 10 -2Q 
for the observation. The burst is shown in Figure 5. 

We note, however, that the field of view contains other 
known magnetars including AX J1845.0—0300 (Torii 
et al. 1998) located 0.38 ° from the center of the point¬ 
ing, and IE 1841—045 (Vasisht & Gotthelf 1997), located 
2.3° from the center of the pointing. As IE 1841—045 
is an active and frequent burster (e.g. Lin et al. 2011; 
An et al. 2015), it is possible that the burst originated 
from this source. While we cannot exclude the possibil- 



Fig. 5.— Burst from the direction of PSR J1846—0258 on MJD 
55070. The time series is binned with 1-ms time resolution and 
covers the full 2-60 keV range of the PCA of both operational 
PCUs. 

ity that this burst originated from PSR J1846—0258, we 
note that there is no change to either the radiative prop¬ 
erties, or timing behavior at these epochs to within our 
measurement uncertainties. 

5. DISCUSSION 

We have presented seven years of post-outburst timing 
of PSR J1846—0258 in which we measure the braking 
index to be n = 2.19 ± 0.03. This is discrepant at the 
14.5<7 level from the pre-outburst braking index of n = 
2.65 ± 0.01 (Livingstone et al. 2006). We note that this 
measurement is made over a comparable span of time 
to that over which the pre-outburst braking index was 
measured. 

Only one other rotation-powered pulsar has 
had a radiative change associated with a glitch: 
PSR J1119—6127. Following a glitch in 2007, the radio 
pulse profile changed from single- to double-peaked. 
This double-peaked profile was only seen once, during 
the first post-glitch observation of the pulsar, and 
had returned to the single-peaked profile by the next 
observation (Weltevrede et al. 2011). It appears that 
PSR J1119—6127 may have undergone a change in 
braking index of similar magnitude following this 
radiatively loud glitch, with a ~ 15% reduction in n 
at the time of the glitch (Antonopoulou et al. 2015). 
However, only formal phase-connected timing errors are 
given for this possible change in the braking index, and 
this method is susceptible to timing noise, (e.g. Hobbs 
et al. 2010; Livingstone et al. 2011). Given this, and the 
large non-white residuals seen after the fitting, the true 
significance of this result is currently unknown. 

It is interesting that the only two nominally rotation- 
powered pulsars which have been observed to have ra¬ 
diatively loud glitches are two of those with the high¬ 
est dipole-inferred magnetic field. In both cases the ob¬ 
served braking indices were consistent with being con¬ 
stant through radiatively quiet glitches and decreased 
following their loud glitches. This decrease in brak¬ 
ing index effectively has the pulsars moving faster to¬ 
wards the magnetar population on the P-P diagram. 
This, together with radiatively loud glitches being a 
defining characteristic of magnetars (e.g. Dib & Kaspi 
2014) is suggestive that the large magnetic field in these 
two seemingly rotation-powered pulsars is responsible for 
their unusual activity. 

There was also a change in v in the high-magnetic- 
field rotation-powered pulsar PSR J1718—3718 follow- 
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ing a large glitch (Manchester & Hobbs 2011). While 
the implied v both before and after this pulsar’s glitch 
gives nonphysical braking indices, n ~ —17(5) and 
n ~ —146(2), the measured i) were consistent over 
~ 3000 days before the glitch, and for the ~ 700 days 
after it. Again, while the implied braking indices seem 
nonphysical, it is interesting that D changed with a glitch 
in yet another high-magnetic-field pulsar. 

A possible change in the braking index was seen in 
the Crab pulsar, where for a ~ 11-yr span the measured 
braking index was ~ 8% lower than the long-term av¬ 
erage braking index. This period of low braking index 
occurred during a period of higher-than-normal glitch 
activity, and Lyne et al. (2015) note that this possible 
change is most likely due to unmodeled glitch parame¬ 
ters. 

One possibility to explain a substantial change in 
a braking index, such as the one we observe in 
PSR J1846—0258, would be contamination due to a long¬ 
term glitch recovery. If this is the case, one would expect 
a bias towards a higher v. and thus a higher n (Lyne et al. 
2000). This is due to the typical glitch behavior of an 
exponentially decaying v, which leads to a decrease in 
the magnitude of the measured v as a function of time, 
and thus to an artificially larger braking index. This is 
the opposite of what we observe. 

There are several theoretical models to explain the 
observation that all measured braking indices are less 
than the canonical n = 3 of a magnetic dipole in a vac¬ 
uum. As yet, the change in braking index observed in 
PSR J1846—0258 is unique - it is larger than ever before 
seen, and appears to be constant following the magnetar- 
like event. Here we will discuss the consequences of a 
changing braking index in the context of these models. 

In particle-wind models (see e.g. Harding et al. 1999; 
Tong et al. 2013), one can explain any braking index 
between n = 1 — 3 by combining spin-down effects from 
both the standard magnetic dipole radiation (n = 3) 
with that of angular momentum loss from an out flowing 
particle wind (n = 1). As shown in Lyne et al. (2015), 
one can express the fraction of spin-down power due to 
a particle wind as: 


This would imply that before 2006, 21±1% of 

PSR J1846—0258’s spin-down was due to a wind, and 
68 ± 4% after 2006. This model predicts a relation be¬ 
tween the braking index and the luminosity of the parti¬ 
cle wind (Harding et al. 1999; Livingstone et al. 2011): 


L p = (3 - n) 2 



6/ 2 c 3 

B 2 R 6 


( 6 ) 


uncertainties, ~10%. 

One can also obtain a braking index different from 3 
by relaxing the assumption of a constant magnetic dipole 
in a vacuum, allowing the dipole to change over time, 
(see e.g. Gunn & Ostriker 1969; Manchester et al. 1985; 
Blandford & Romani 1988). This is expressed in a con¬ 
venient form by Lyne et al. (2015): 


v 

flobs = Tldip T V 
V 



a 

tana 



(7) 


To explain a braking index lower than ridip of 3, either 
the moment of inertia I is decreasing, or either the mass 
M or the angle of mis-alignment between the spin and 
magnetic axis a are increasing. Furthermore, this implies 
that at the epoch of the magnetar outburst, the fractional 
rate of change of the magnitude of either I/I, d/tana, 
or M/M increased by a factor of 2.3±0.2. It does not 
seem physically plausible to have so large a change in 
either I/I nor M/M, especially given the lack of change 
of the flux of the pulsar wind nebula. Such a change in 
either d or a also seems improbable, given the lack of 
any detected change in the pulse profile (see §4.2). 

One could also change the braking index by altering 
the geometry of the magnetosphere (see e.g. Thompson 
et al. 2002; Contopoulos & Spitkovsky 2006). In the 
twisted neutron-star magnetosphere model of Thomp¬ 
son et al. (2002), the braking index of a pulsar is given 
by n = 2p + 1, where p is radial index. The observed 
change in braking index in this model implies that for 
PSR J1846—0258, the “twist” between the north and 
south hemispheres increased by ~ one radian at the time 
of the outburst, which should lead to a corresponding in¬ 
crease in the X-ray luminosity of ~ 50%. This is not 
seen. Additionally, in any magnetospheric origin for a 
change in braking index, one would need to modify the 
magnetosphere in such a way as to maintain a constant 
pulse profile over the magnetar-like event, which seems 
challenging. 

Beloborodov (2009) has a modified version of this 
model in which instead of a global twist in the mag¬ 
netosphere, the twist is concentrated into a localized re¬ 
gion known as a “j-bundle.” This j-bundle will increase 
the dipole moment of the neutron star, leading to an 
increased spin-down rate. As the j-bundle shrinks, the 
effective dipole moment should decrease with time, lead¬ 
ing to a positive contribution to u and thus the brak¬ 
ing index. While this model can be used to explain the 
glitch behavior of PSR J1846—0258 associated with the 
magnetar-like event (Livingstone et al. 2010), it does not 
immediately explain our observed long-term decrease in 
braking index. 


where / is the moment of inertia, B the magnetic field, 
and R the pulsar’s radius. Assuming neither the mag¬ 
netic field nor moment of inertia changed substantially, 
the luminosity of the pulsar wind nebula might have been 
expected to increase by a factor of approximately 5. Such 
a significant flux change was ruled out by deep Chandra 
observations by Livingstone et al. (2011), as well as by 
the consistency of the flux during the Swift campaign 
with the pre-outburst flux reported by Kumar & Safi- 
Harb (2008) to within the telescopes’ cross calibration 


6. CONCLUSIONS 

The observed braking index of PSR J1846—0258 has 
significantly changed following its period of magnetar- 
like behavior. This long-term change in n is, to within 
measurement errors, unaccompanied by any correspond¬ 
ing long-lived change in the flux of the source, or any 
change in its pulse profile. This is in contrast to most 
of the models discussed above where a correlated change 
in the X-ray luminosity is expected for both wind-based 
models (Harding et al. 1999) and global magnetospheric 
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twist based models (Thompson et al. 2002). As well, 
models which modify the assumptions of a constant mag¬ 
netic dipole require far too high a change in / or M to 
be physically plausible, or a change in a or d which seem 
unlikely given the stable pulse profile. 

The most plausible explanation for a changed braking 
index appears to be due to some form of change in mag- 
netospheric configuration, but this change is constrained 
by our observations to be unaccompanied by any large- 
scale change in flux, spectrum or pulse profile. One pos¬ 
sible way to probe the magnetosphere of pulsars, and 
therefore test this hypothesis, is by means of X-ray po- 
larimetry. Measurements of polarization fractions and 
angles are very sensitive to viewing geometries, as well 
as twists in the magnetosphere (e.g. van Adelsberg & Lai 
2006; Taverna et al. 2014). 
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